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Correlations of eigenfunctions, {\ipk(ri)\ 2 \ipi(r2)\ 2 } , in a 
disordered system are investigated. We derive general for- 
mulae expressing these correlation functions in terms of the 
supcrmatrix sigma-model. In particular case of weak localiza- 
tion regime we find that the correlations of the same eigen- 
function are proportional to g~ for large distances, while the 
correlations of two different eigenfunctions cross over from 
<? _1 behavior for r\ = T2 to g~ 2 one for [n — ra| 2> /, with g 
and / being the dimensionless conductance and the mean free 
path, respectively. 

PACS numbers: 05.45+b, 72.15.-v, 73.20.Dx 



Statistics of eigenfunction fluctuations in disordered 
and chaotic systems have attracted a research interest re- 
cently. The fluctuations of eigenfunction amplitudes de- 
termine statistical properties of conductance peaks and 
level width in quantum dots in the Coulomb blockade 
regime 0-^], and can be directly measured in the mi- 
crowave cavity experiments . On the theoretical side, 
the recent progress is based on application of the super- 
symmetry method to the problem of eigenfunction statis- 
tics [|]]^] . ft was found that the distribution of eigenfunc- 
tion amplitudes is in the leading (zero-mode) approxima- 
tion correctly described by formulae of the random ma- 
trix theory (RMT). Deviations from the RMT predictions 
were studied in Refs. |^-|To[]. Correlations of amplitudes 
of an eigenfunction in two different spatial points were 
considered in Ref. [jll| on the level of zero-mode approx- 
imation; the latter was shown fl2] | to be equivalent to 
the RMT-like assumption of the Gaussian fluctuations of 
wavefunctions. 

All results mentioned above concern fluctuations of the 
same eigenfunction. In the present Letter we study for 
the first time correlations of amplitudes of two different 
eigenfunctions. We derive general expressions in terms of 
the supermatrix sigma model, valid for arbitrary diffusive 
(or classically chaotic) system, and then apply them to 
the weak localization regime. 

In order to evaluate the correlations of the wavefunc- 
tions we use the technique similar to that of Ref. [p"3| . 
Namely, we consider a quantity 

A(r u r 2 ,uj) (1) 
= ( E \Mri)Mr2)\ 2 S(e - e k )5(e + 



Here the angular brackets denote the impurity average. 
We have introduced the eigenstates ip k (r) and eigenval- 
ues e k of the Hamiltonian H = Hq + U (r) in a particular 
disorder configuration U(r), Hq being the Hamiltonian 
of the free particle. Further we denote 

a(n,r 2 ,e k ) = (\ip k (ri)i() k {r 2 )\ 2 ) 
P(r ll r 2l e kl e l ) = (\Mri)Mr2)\ 2 ), M *> (2) 

where the averaging is carried over all eigenstates ip k i ipi 
with given values of the energies e k , £i- Since we are in- 
terested in a relatively narrow window in the spectrum, 
lu <C €f, these quantities are translationally invariant: 



a(r 1 ,r 2 ,e k ) = a(ri,r 2 ), 
P{ri,r 2 ,e k ,ei) = f3(r 1 ,r 2 ,e k - £;)■ 



(3) 



Due to the presence of delta- functions in Eq. (|l|) , we can 
rewrite the latter exactly as 



A(r 1 ,r 2 ,iv) = a{rx,r 2 )A 1 S(u) 

+ /3(r lir2 ,u;)A- 2 R(Lu), 



(4) 



where A is the mean level spacing, A = (vV)^ 1 , with 
V and v being the system volume and the density of 
states, respectively. Here we have introduced the two- 
level correlation function, 



R(oj) 



]S(e- e k )S(e + uj - e t ) 



(5) 



We would like to stress that all transformations up to 
now are completely rigorous. In particular, we do not 
assume any decoupling of eigenfunction and eigenvalue 
correlations. 

On the other hand, the quantity (|l|) can be writ- 
ten in terms of the Green's functions in the coordinate- 
frequency representation, 

A(ri,r 2> w) = (2TT 2 )- 1 
xRe { (G R (r 1 ,r 1 ,e)G A (r 2 , r 2 ,e + u)) 
- (G R (r 1 ,r 1 ,e))(G A (r 2 ,r 2 ,e + cj))}. (6) 

The expression (||) can be directly calculated with the 
use of the supersymmetry technique Jl4|-^6[. We con- 
centrate in the sequel on the case of broken time-reversal 
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symmetry (unitary ensemble) ; generalization to the other 
ensembles is straightforward. After the standard manip- 
ulations we get 

A(ri,r 2 ,w) = -(2 7 r 2 )- 1 Re{< 5h 1 , 1 (r 1 ,r 1 ). 9b 2 b 2 (r 2 ,r 2 ) 
+ 9bb( r u r 2)gll(r2,ri)) F - (gll{r 1 ,r 1 )) F 
x{g 22 (r 2 ,r 2 )) F }. (7) 

Here (. . .) F denotes the averaging with the action of the 
supermatrix sigma-model i^[<5]: 

(...)r = ! DQ(...)exp(-F[Q]), 

m--*j ^^mm, m 

where D is the diffusion coefficient, Q = T _1 AT is a 
4x4 supermatrix, A = diag(l, 1,-1,-1), and T belongs 
to the supercoset space U(l, l|2)/f7(l|l) x U(l\l). The 
symbol Str denotes the supertrace defined as Str_B = 



nil nil 
B bb ~ B ff 



Bbb ~ B 22 . The upper matrix indices corre- 



spond to the retarded-advanced decomposition, while the 
lower indices denote the boson-fermion one. The Green's 
function g in Eq. (Q) is the solution to the matrix equa- 
tion: 



-i(e + | - H ) - % -{to + »0)A + Q/2t 
5(n - r 2 ). 



g(ri,r 2 ) 



(9) 



Expressing these functions through the matrices Q and 
taking into account Eq. (Q), we arrive at the following 
equation valid in arbitrary diffusive system: 



2tt 



6(u) 



R(u) 



A ' A 2 

-(^) 2 Re(Qi b 1 (r 1 )Q 22 (r 2 )) F 
[Im G R ( ri - r 2 )] 2 Re(Qi 2 (r 1 )g 2 b 1 (r 1 )) 



(to,) 2 , (10) 



with G R being the impurity averaged retarded Green's 
function. In particular, in the case of 2D and 3D system, 
is given by 



G R (r) 



r/2 



tt/2 



d6exp[(ip F r -r/2l) cos 6], 2D 



—irv(p F r) 1 e~xjp[ippr — r/2l], 



3D 



where I is the mean free path. The key point making 
the further progress possible is that the term containing 
the single eigenfunction correlations in lhs of Eq.(|lO[) is 
proportional to S(oj), whereas the one depending on the 
correlations of two different eigenfunctions is regular at 
lu = 0. Thus, separation of the expression in the rhs 
of Eq. jicj ) into the singular (proportional to 6(ui)) and 
regular parts allows one to obtain the quantities a(r\, r 2 ) 
and P(ri,r 2 , a). 



Now we turn to the case of a metallic system in the 
weak localization regime. The corresponding small pa- 
rameter is given by Eq. (|l^) below. For further purposes, 
we introduce the functions 

/ 1 (r 1 ,r 2 )=n 2 (r 1 ,r 2 ), 

/ 2 (ri,r 2 ) - (2V)- 1 J dr [n 2 (r, n) + n 2 (r, r 2 )] , 
h = V- 2 J drdr'n 2 (r,r'), 



h(ri,r 2 ) = V- 1 J drU(r, ri )U(r,r 2 ). 



(11) 



Here the diffusion propagator n is the solution to the 
diffusion equation 

-£>V 2 n(n,r 2 ) = (Tr^)-M(ri -r 2 ) (12) 

with appropriate boundary conditions. We obtain: 

n(n, r 2 ) = (ttvV)- 1 Y^W)~ l Mri)<f> q fo), (I 3 ) 
q 

with (j) q being the eigenfunction of the diffusion operator 
corresponding to the eigenvalue Dq 2 , q ^ 0. The level 
correlation function has the form fllTI] 



R(lu) = 1 - s~ 2 sin 2 s + f 3 sin 2 s + Q{g~ 



(14) 



where a dimensionless parameter s = 7rw/A is intro- 
duced. The first two terms in Eq. (|l4|) are given by 
RMT, while the third one is the correction of order g~ 2 
due to the diffusion modes. Here g = 2irE c /A is the 
dimensionless conductance, with E c being the Thouless 
energy. 

The sigma-model correlation functions 
{QH(ri)Qll(r 2 )) F and {QU^QU^)) F can be cal- 
culated for relatively low frequencies lu <C E c with the 
use of a general method developed in Refs. which 
allows one to take into account spatial variations of the 
field Q. The results are obtained in form of expansions 
in g~ 
9 



(Qll(ri)QH(r 2 )) F 

= -1 



exp(is)sins 

2« — ; 2t- 



i0 



n(n,r 2 ). (15) 



The first two terms in Eq. ( Jl 5|) represent the result of the 
so-called zero-mode approximation |14| ] , which takes into 
account only the spatially constant configurations of the 
field Q(r), so that the functional integral over DQ(r) 
is reduced to an integral over a single matrix Q. The 
last term is the correction of order g^ 1 . An analogous 
calculation for the second correlator yields |Tq] : 



{Qlt{r l )QH{r 2 )) F = -2 
.exp(is) sins 



1 



(s + toy 



s + iO 

n(n,r 2 ) 



(16) 



2 



Now, separating regular and singular parts in rhs of 
Eq. (|1C|), we obtain the following result for the autocor- 
relations of the same eigenfunction: 



V 2 (\Mri)Mr2)\ 2 )-l 

= k d (r)[l + U(r 1 ,r 1 )]+U(r 1 ,r 2 ), 



(17) 



and for the correlation of amplitudes of two different 
eigenfunctions 

V 2 (\Mri)Mr2)\ 2 )-l = k d (r)U(r 1 ,r 1 ), k^l (18) 

Here r — \r\ — r 2 \, and the function k d (r) is defined as 

k d (r) = {irv)- 2 [lmG R (r)] 2 = 

1, 2D, p F r<l 

exp(— r/l) { Sinppr)^ 1 cos 2 pft, 2D, p F r > 1 
(pFf) ~ 2 sin 2 pFf, 3D 



Note that the result Q17|) for n = r 2 is the inverse partici- 
pation ratio previously obtained in Ref . || , while that for 
arbitrary spatial separation was found in the zero-mode 
approximation (g = oo) in Ref. |ll[| . 

Eq. ( |l8| ) shows that the correlations between different 
eigenfunctions are relatively small in the weak disorder 
regime. Indeed, they are proportional to the small pa- 
rameter II(r,r), which is equal in the cases of 2D and 
3D geometry to (L is the size of the system) 



II(r,r) 



(irg^lnL/l, 
- X L/l, 
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2D 
3D 



(19) 



In particular, for r\ = r 2 we have 

V 2 (\Mr)ipir)\ 2 ) ~ 1 = hi + (1 + 5 k i)IL(r,r). (20) 

The correlations are enhanced by disorder; when the sys- 
tem approaches the mobility edge in 3D or the strong 
localization regime in 2D, the relative magnitude of cor- 
relations, II(r, r) becomes the quantity of order of unity. 

An inspection of Eqs. (|l7|), (18]) shows that while the 
correlations of amplitude of the same wavefunction sur- 
vives for the large separation between the points, r S> I, 
and is proportional to g~ 1 , the correlations of two differ- 
ent wavefunctions decay exponentially for the distances 
larger than the mean free path I. This is, however, an 
artifact of the g^ 1 approximation, and the investigation 
of the corresponding tails requires the extension of the 
above calculation to the terms proportional to g~ 2 . We 
find that the correlator (Q\IQ 22 )f gets the following cor- 
rection: 



HQUQIDf = -fi + 2/4 
.exp(2is 



exp(2is)/ 3 



-2i- 

s + iO 

exp(2is) - 



- 1 



2(s + i0) 2 



h) 

(h- 4/2 + 3/3-4/4). 



(21) 



Consequently, we obtain the following results for the cor- 
relations of different eigenfunctions at r > I: 



V 2 (\Mri)Mr2)\ 2 )-l 



h-h- 2/4 



1 ( sin 2 s x 



2 v 
2sin(2s) 



(/ 2 - h) 



sin 2 s 



(fi - 4/ 2 + 3/3 - 2/4 



k ^ I (22) 



As it should be expected, the integrals over the coordi- 
nate of this correlation function are equal to zero. This 
property is just the normalization condition and should 
hold in arbitrary order of expansion in g . Note also 
that the correlation of different eigenfunctions becomes 
now frequency dependent in a non-trivial way, even in 
the regime u> -C E c considered. 

The quantities f 2 , fy, and fa are proportional to g~ 2 , 
with some numerical prefactors Jl7| . On the other hand, 
/1 depends essentially on the distance r = |r*i — r 2 \. In 
particular, for I C r < L we find 



/ 1 (r 1 ,r 2 )=n 2 (r 1 ,r 2 ) 



(4n 2 vDr) 2 ' 



2D 
3D 



Thus, for I < r <C L, the contributions proportional to 
/1 dominate in Eq.(|22|), and we get 

V^(ri)^(r-2)| 2 )-1 = V(r 1; r 2 ) , k ± I. (23) 

Another correlation function, generally used for the 
calculation of the linear response of the system, 

-f(ri,r 2 ,u>) = {ip* k {ri)tl)i{ri)^k{r 2 )^i {r 2 )) , k^l, 

can be calculated in a similar way (cf. Starting 
from the quantity 

B(ri,r 2 ,uj) = 

( ^i{rx)^l{r x )^ k {r 2 )^{r 2 )8{e - e k )5(e + lj - q) ) 



k,l 



^2tpl(n)ip k (r 2 )S(e - e fc ) ) 
\ k I 

^i(ri)^(r 2 )S(e + uj - a) ) , 



(24) 



and repeating the derivation that led us to Eq. ([To|), we 
get another identity: 



2tt 



a(rx,r 2 ) 



8(u) 



7(ri,r 2 ,w) 



R(uj) 



A ~ v ~' ' A 2 
-M 2 Re{(g b 12 (r 1 )g 2 b 1 (r 2 )) i , 
Ur^QlK^QlK^F ~k d {r)}, 



(25) 
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Taking into account Eqs. ( |l5| ) and (jiff), and separating 
the rhs into the regular and singular parts, we recover 
Eq. (|l^) and obtain 



(Y. M. B.) and SFB195 der Deutschen Forschungsge- 
meinschaft (A. D. M.). 



V*(^(riM(ri)^(raM*(ra)> 
= k d {r)+H(r 1 ,r 2 ), k^l. 



(26) 



As was mentioned, the above derivation is valid for 
lu <C E c . In order to obtain the results in the range 
lu > E c one can calculate the sigma-model correlation 
functions entering Eqs. (fTo|), j25| ) by means of the per- 
turbation theory |p|] . We find then for k I 

V 2 (\Mri)Mr2)\ 2 ) = 1 + Re{k d (r)n u (r 1 ,r 2 ) 

+ \ nS(ri,ra)-i j drdr'Ul(r,r') }, (27) 

V 2 m(ri)MrMr2)^t(r 2 )) = k d (r) + Ren^n, r 2 ), 

where n^(ri,r 2 ) is the finite-frequency diffusion propa- 
gator 



4>q{ri)4>g(r2) 
Dq 2 — iu> 



(28) 



and the summation in Eq. ( p8| ) now includes q — 0. 

The above derivation was performed for the case of 
a disordered sample. However, one can repeat it for 
a classically chaotic ballistic system, using the recent 
derivation of the supersymmetric sigma-model for this 
case Jl^]. In the end one gets the analogous results, but 
with the diffusion operator being replaced by the Perron- 
Frobenius one for the given chaotic system. To analyze 
the results quantitatively in this case, one needs then an 
information about the eigenfunctions and eigenvalues of 
the Perron- Frobenius operator. This question deserves 
further investigation and goes beyond the scope of this 
Letter. 

Finally, we discuss the relation between our results and 
the recent experiment ||, where strong correlations in 
amplitudes of neighboring conductance peaks of a quan- 
tum dot were observed. Within the one-electron picture, 
this seems to contradict our result ( |l8| ) of weakness of 
correlations of different eigenfunctions. In principle, one 
could imagine that the dot was far from the universal 
(RMT) regime, so that the parameter H(r,r) determin- 
ing the magnitude of correlations was not small. This 
would be however in contradiction with the fact that the 
total distribution of peak heights in Ref. § was well de- 
scribed by the RMT formulae jjj, since the corrections 
to the distribution of |'0 2 (r)| (and consequently to that 
of peak heights) are proportional to the same parame- 
ter II(r,r) g|. Another possibility is that the electron- 
electron interaction effects lead to some modification of 
eigenfunction correlations. This problem remains to be 
studied in future. 

We are grateful to Charlie Marcus for a discussion 
which stimulated this research. The work was sup- 
ported by the Alexander von Humboldt Foundation 
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